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Multiplicative processes and multifractals have earned increased popularity in applica- 
tions ranging from hydrodynamic turbulence to computer network traffic, from image 
processing to economics. We analyse the multifractality of the recently proposed point 
process models generating the signals exhibiting noise. The models may be used 

for modeling and analysis of stochastic processes in different systems. We show that 
the multiplicative point process models generate multifractal signals, in contrast to the 
formally constructed signals with l//' 3 noise and signals consisting of sum of the uncor- 
rected components with a wide-range distribution of the relaxation times. 



1. Introduction 

Multifractal models are used to account for scale invariance properties of various 
objects in different domains ranging from the energy dissipation in turbulent flowJU 
to financial datcPl. Healthy human heartbeat intervals exhibit multifractal proper- 
ties rather than being fractal for a life-threatening condition, known as congestive 
heart failure^. Cerebral blood flow in healthy humans is also multifractal^. 

Scaling behavior has become a welcome careful description of complexity in 
many fields including natural phenomena, human heart rhythm in biology, spatial 
repartition of faults in geology, as well as human activities such as traffic in computer 
networks and financial markets. The multifractal formalism has received much 
attention as one of the most popular frameworks to describe and analyse signals 
and processes that exhibit scaling properties, covering and connecting both the local 
scaling and the global one in terms of sample moments. 

The purpose of this paper is to analyse the multifractality of signals exhibiting 
1/fP noise generated by differen t tec hniques and, especially, of the point processes 
with 1/ fP power spectral densit jEG2. 

First of all, however, we will analyse the multifractality of the signal constructed 

[71 

by the inverse fast Fourier transform . Using this method we can generate signals 
with any desirable slope (3 of the power spectral density S(f) ~ f//' 3 - 

We calculate a generalized qth order height-height correlation function (GHCF) 
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F q (t) defined aM 

F q (t) = (\i(t'+t)-i(t')n 1/q , a) 

where the angular brackets denote the time average. The GHCF F q (t) characterizes 
the correlation properties of the signal I(t), and for a multiaffinc signal a power-law 
behavior like 



F q (t)~t H < 



(2) 



is expected. Here H q is the generalized gth order Hurst exponent. If H q is inde- 
pendent on q, a single scaling exponent H q is involved and the signal I{t) is said to 
be monofractaM If H q depends on g, the signal is considered to be multifractal. 
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Figure 1. a) Power spectral density and b) the generalized Hurst exponents H q versus 1/q in the 
scaling regime 1 < t < 1000 for the slopes /3 = 1 (open circles) and /3 = 2 (open squires). The 
signal of 10 e points was generated and averaged over 10 realizations, c) and d) show GHCF F q (t) 
versus time t for the same parameters /3 = 1 and [3 = 2, respectively. 



In figure [T] a) we present a power spectral densities with the different slopes (3 
and in figure[T]b) we show the Hurst exponents, calculated from GHCF using linear 
regression dependence on 1/q of the signals formally constructed by the inverse 
Fourier transform. In figure Q]c) and d) corresponding GHCF F q (t) versus time t 
are shown. We see that Hurst exponent H q does not depend on q, which indicates 
that the signal is monofractal. 
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2. Stochastic multiplicative point process 

In many cases the intensity of some signals or currents can be represented by a 
sequence of random (however, as a rule, mutually correlated) pulses or elementary 
events Ak(t — tk), 

i(t) = J2Mt-tk). (3) 

Here the function Ak(4>) represents the shape of the k pulse making an influence 
on the signal I(t) in the region of the transit time t k - We will be interested in 
the processes with the power-law distribution of the power spectral density at low 
frequencies. It is easy to show that the shapes of the pulses mainly influence the 
high frequency, / > 1/At p , with At p being the characteristic pulse length, power 
spectral density, while fluctuations of the pulse amplitudes result, as a rule, in the 
white or the Lorentzian but not 1// noised Therefore, we restrict our analysis to 
noise due to correlations between the transit times tk- In such approach we can 
replace the function Ak(t — tk) by the Dirac delta function and then express the 
signal as 

J(t) = o £ *(*-**)> (4) 
fc 

with a being an average contribution to the signal of one pulse. This mo deP also 
corresponds to the flow of identical objects: electrons, photons, cars, and so on, and 
is called the point process model. Point processes arise in different fields, such as 
physics, economics, cosmology, ecology, neurology, seis molog y, traffic flow, signaling 
and telecom networks, and the Internet (see e.g., paper a I land references herein). 

The power spectrum of the point process signal is described completely by the 
set of the interevent intervals Tk — ifc+i — tk- Moreover, the low frequency noise is 
defined by the statistical properties of the signal at a large-time-scale, i.e., by the 
fluctuations of the time difference 

k+q-l 

A(k; q) = t k+q - t k = ^ Ti ( 5 ) 

i—k 

at large q, determined by the slow dynamics of the average interpulse time 
fk(q) — A(k;q)/q between the occurrence of pulses k and k + q. Quite generally 
the dependence of the average interevent time ffc may be described by the general 
Langevin equation. The Langevin equation may be written down in the actual time 
t or, equivalently, in the space of the occurrence numbers k with the drift coefficient 
h(fk) and a multiplicative noise g(fk)^(k), 

^ = h(f k )+g(f k )m- (6) 

Here we interpret fc as a continuous variable while the white Gaussian noise £(fc) 
satisfies the standard condition 

m)ak'))=5(k-k') (7) 
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with the brackets (. . .) denoting the averaging over the realizations of the process. 
We understand the equation J6| in Ito interpretation. 

Transition from the occurrence nu mbe rs k to the actual time t in Eq. J6]) may 
be fulfilled using the relation dt = ffcdAPl 

The particular sequence of the interevent times Tk may be superimposed by 
some additional noise or stochasticity, e.g., Tk may be determined by the Poisson 
distribution 

P(r k ) = ^e- T «l"« (8) 

Tk 

with the slowly, according to Eq. ^ , changeable average interevent time f ^ . Such 
additional stochasticity do not influence the long-range statistical properties and 
the low frequency spectra of the process. Therefore, further we will restrict the 
analysis to the processes generated by Eq. ^ and will identify Tk with ffc. 



2.1. Power spectral density 

The point process is entirely defined by the occurrence times tk ■ The power spectral 
density of the point process Q may be expressed as 

2 



S(f) = lim ( - 



I{t)e- t27Tft dt 



\ k q—k m i n —k i 



(9) 

where ti and tf are initial and final observation times, T = tf — t% 3> lu^ 1 is the 
whole observation time and u> = 2irf. Here fc m i n and fc max are minimal and maximal 
values of index k in the interval of observation T and the brackets (. . .) denote the 
averaging over realizations of the process. 

For the interpulse intervals described by the Langevin equation ([6]) we use a 
perturbative solution in the vicinity of Tk ■ After replacing the averaging over k by 
the averaging over the distribution Pk{Tk) of the interpulse times Tfe, we have the 
power spectrum^ 

oo 

S(f) = 2l 2 ^ 7 / Pk(T k ) Re[e-^-T) erfc v^]^^, (10) 



where I and f are the averages of the signal and the interpulse times, respectively, 
and x = 7r/r|//i(rfc). 

The replacement of the averaging over k and over realizations of the process by 
the averaging over the distribution of the interpulse times Tk, Pk (Tk), is possible 
when the process is ergodic. Ergodicity is usually a common feature of the stationary 
process described by the general Langevin equation. 

According to Eq. (jTU]) the small interpulse times and the clustering of the pulses 
make the greatest contribution to 1/ f& noise. The power-law spectral density is 
very often related with the power-law behavior of other characteristics of the signal, 
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such as autocorrelation function, probability densities and other statistics, and with 
the fractality of the signals, in generaff^l. Therefore, we investigate the power-law 
dependences of the drift coefficient and of the distribution density on the time Tfc 
in some interval of the small interpulse times, i.e., 



h{r k ) = Pk(T k ) = Ct£, T min < T fc < T max , (11) 

where the coefficient 7 represents the rate of the signal's nonlinear relaxation and 
C has to be defined from the normalization. 

The simplest and the well-known process generating the power-law probability 
distribution function for r k is a multiplicative stochastic process with g(r k ) = art 
and S = 2fj, — 1, written a d^^ 

Tk+l =T k + JT^' 1 + <7T^e k . (12) 

Here 7 represents the nonlinear relaxation of the signal, while T k fluctuates due to 
the perturbation by normally distributed uncorrelated random variables e k with a 
zero expectation and unit variance and a is a standard deviation of the white noise. 

Eq. (|T2"|) is the difference (discrete) version of the differential equation ((6|) . On 
the other hand, it is the generalization of the simple autoregressive model of 1// 
noise^ (see also Eq. ([15])) and represents quite general evolution of the interevent 
time with the nonlinear drift h(r k ) = 1 and the multiplicative noise ar^e k , 

resulting in the l//' 3 noise and power-law distribution (|lljl of the interevent time 
Tfc with the exponent a = 27/cr 2 — 2/i. Indeed, the power spectrum for the process 
(HU), when 7/(7^) « / « 7 /(t^), 

(2 + a )(/3-l)a 2 r(/3-l/2) / 2 n/?-i 1 
* UJ ^a{rl + « - r 2 | n «) sin(7r/3/2) \J ff ( 

where 

a=^-2/i, 0=1+ —2L-, l</3<2. (14) 
a A 6 — 2/i 2 

For fj, = 1 we have a completely multiplicative point process when the stochastic 
change of the interpulse time is proportional to itself. Another case of interest 
concerns /1 = 1/2, then we have the Brownian motion of the interevent time with 
the linear relaxation of the signal I ~ a/r. 

Figure [2] represents the spectral densities (j9]) with different slopes /3 of the signals 
generated numerically according to Eqs. and (Tj"2"|) for different parameters of the 
model. We see that the simple iterative equation (fT2|) with the multiplicative noise 
produces the signals with the power spectral density of different slopes, depending 
on the parameters of the model. The agreement of the numerical results with the 
approximate theory is quite good. 
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Figure 2. Power spectral density (J9j vs frequency of the signal generated by Eqs. Q and (1 12t 
with the parameters a) fi = 0.5, cr = 0.02 and different relaxations of the signal 7 = 0.0001 (open 
circles), 0.0002 (open squares) and 0.0003 (open triangles); and b) /i = 1, a = 0.1 and 7 = 0.008 
(open circles), 0.01 (open squares) and 0.012 (open triangles). We restrict the diffusion of the 
interevent time in the interval r m ; n = 10 -6 , T max = 1 with the reflective boundary condition at 
-r min and transition to the white + crej, , for Tfc > T max and 100 realizations 

with 10 6 tfc points each were used. The solid lines represent the analytical results according to 
Eq. | Tl3t . 



2.2. Multifractal point processes 

The multifractal formalism has received much attention recently as one of the most 
popular frameworks to describe and analyse signals and processes that exhibit the 
scaling properties. 

Fractality of the point process can be investigated by transition from the point 
process to the stochastic signal I(t), using the rectangular constant area pulses, 
instead of the Dirac delta functions. The stochastic signal will have the same fractal 
properties as the origin point process. 




L0" 1 10" 10 1 10 2 10 3 10 4 10 5 0.5 1 1.5 



t l/q 

Figure 3. a) Generalized height-height correlation function F q (t) versus time t for the signal JJ) 
and JT2} with fi = 0.5, f = 1, a = 0.02, 7 = 0.0002, and g = 0.5, 1, 1.5, . . . , 5 from bottom to 
top. b) The generalized Hurst exponents H q versus l/q in the scaling regime 1 < t < 1000 for 
the f = 1 and fj. = 0.5, a = 0.02, 7 = 0.0002 (open circles); n = 0.5, a = 0.02, 7 = 0.0003 (open 
squares); fi = 1, a = 0.1, 7 = 0.008 (open triangles). 
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In figure [3] a) we present the GHCF as a function of the time interval t, and in 
figure [3] b) we show the Hurst exponents calculated from GHCF using the linear 
regression dependence on 1/q for different parameters jjL, a, and 7. We observe 
the clear multifractal behavior since the slopes of the log-log plot of GHCF are 
depending on q. 

Another interesting case is a sequence of transit times with random increments 
of the time intervals between pulses, Tk = Tk-i + crsk- It is natural to restrict in 
some way the infinite Brownian increase or decrease of the interpulse times Tk, e.g., 
by the introduction of the relaxation to the average interpulse time f rate 7. So, 
we have an additive point process 

Tk = Tfc_i - 7(r fc _i - f) + ae k . (15) 

This model generates the process with 1// noisc^, and may be useful for modeling 
and analysing different systems (see references in paper^). Introduction of the 
reflective boundary condition at T m ; n > avoids the formation of clusters and leads 
to l// 2 noise. 
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Figure 4. a) Power spectral density and b) the generalized Hurst exponents of the additive point 
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process ||15|I with the reflecting boundary condition at T mln = 0.1 in the scaling 
t < 1000. The signal of 10 6 points was generated with the parameters f = 1, a 
7 = 0.000001. 



regime 1 < 
= 0.001, and 



In figure [4] a) we present a power spectral density of the additive point process 
(|15p with the reflecting boundary condition and in figure [4] b) we show the gener- 
alized Hurst exponents. We observe the multifractal behavior of the additive point 
process. 



2.3. Monofractality of the white and the Gaussian noises 

It is well-knowii^l that the point processes with the interevent time Tk distributed 
according to Poisson distribution 



P(Tk) 



l e -r k /f 



(16) 
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Figure 5. a) Generalized height-height correlation function F q (t) versus time t, and b) the gen- 
eralized Hurst exponents H q versus 1/q of the point process with interevent time distributed 
according to Poisson distribution l|16j l in the scaling regime 1< t < 1000. The signal of 10 6 points 
with the parameter f = 1 has been generated. 



generate the white noise, S(f) = const. 

In figure [5] a) we present the GHCF as a function of the time interval t and in 
figure [5] b) we show the Hurst exponents for the white noise. The Hurst exponents 
of the white noise are equal to zero. This demonstrates the absence of the scaling 
and that there is no correlation in time. 

Noise with the power-law spectral density 1// is often modeled as the sum of 
the Lorentzian spectra with the appropriate weights of a wide range distribution of 
the relaxation times r rel . The signal may be expressed as a sum of N uncorrelated 
components, 



N Im. 

Z(t) =£/«(*) = / 
1=1 J. 



7(t, 7)5(7)^7, 



(17) 



where 5(7) is the distribution of the relaxation rates 7 
nent J; satisfies the stochastic differential equation 



1 /t , and every compo- 



(18) 



Here Ii is the average value of the signal component J; , (t) is the ^-correlated white 
noise, (£z(t)£z' (t')) — 8i.v5(t~t'), and 01 is the intensity (standard deviation) of the 
white noise. The steady-state solution of the stationary Fokker-Planck equation 
corresponding to Eq. (|18|) for each component Ii and the resulting signal (fl7|) yields 
the Gaussian distribution densities, however, the power spectrum may be of the 
power-law form when a 2 (7)3(7) is constant or a power-law function^. 

In figure [6] a) we present a power spectral density of the sum of the signals with 
a wide range distribution of the relaxation times r rel and in figure [6] b) we show the 
Hurst exponents, calculated from GHCF using linear regression dependence on 1/q. 
In the figures we observe 1// behavior of the signal noise and clearly see that Hurst 
exponent H q does not depend on q, which shows that the signal (fT7|) is monofractal. 
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Figure 6. a) Power spectral density and b) the generalized Hurst exponents H q versus 1/q in the 
scaling regime 1 < t < 1000. The signal of 10 6 points was generated from 10 components with the 
parameters I = 20, a 2 (7)9(7) = 10, and uniform distribution of lg7 with 7 values in the interval 
10~ 4 - 1. 

3. Conclusions 

The multiplicative (fl"!?|) and additive (JTSJ) stochastic point processes may generate 
time series exhibiting the power spectral density S(f) ~ l//' 3 and show clear multi- 
fractal behavior, however, the formally constructed by the inverse Fourier transform 
signal, Poisson white noise (fl"6|) and Gaussian, Eqs. (fT7|) and (fl~8|) . signal with the 
power spectral density S(f)~ 1// are monofractal. 

Therefore, the proposeoP and generalized^ point process models of 1/ ' f 13 noise 
may be used for modeling of stochastic multifractal processes in different systems. 
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